DELAY EQUATIONS DRIVEN BY ROUGH PATHS 



A. NEUENKIRCH, I. NOURDIN AND S. TINDEL 



o 
o 

(N 
> 

o 



Abstract. In this article, we illustrate the flexibility of the algebraic integration formal- 
ism introduced in M. Gubinelli (2004), Controlling Rough Paths, J. Funct. Anal. 216, 
86- 140, by establishing an existence and uniqueness result for delay equations driven by 
rough paths. We then apply our results to the case where the driving path is a fractional 
Brownian motion with Hurst parameter H > \. 



1. Introduction 



In the last years, great efforts have been made to develop a stochastic calculus for 
fractional Brownian motion. The first results gave a rigorous theory for the stochastic 
integration with respect to fractional Brownian motion and established a corresponding 
Ito formula, see e.g. [H El El El [IE]- Thereafter, stochastic differential equations driven 
by fractional Brownian motion have been considered. Here different approaches can be 
P_i ! used depending on the dimension of the equation and the Hurst parameter of the driving 

fractional Brownian motion. In the one-dimensional case [17], existence and uniqueness 
of the solution can be derived by a regularization procedure introduced in [21]. The case 
of a multi-dimensional driving fractional Brownian motion can be treated by means of 
fractional calculus tools, see e.g. [121 [22] or by means of the Young integral [T3|, when 
the Hurst coefficient satisfies H > ^. However, only the rough paths theory [THl and 
^ ' its application to fractional Brownian motion [5J allow to solve fractional SDEs in any 

^ . dimension for a Hurst parameter H > ^. The original rough paths theory developed 

^ I by T. Lyons relies on deeply involved algebraical and analytical tools. Therefore some 

CN I alternative methods [HI E] have been developed recently, trying to catch the essential 

results of [12] with less theoretical apparatus. 



. Since it is based on some rather simple algebraic considerations and an extension of 

^ ' Young's integral, the method given in [0], which we call algebraic integration in the sequel, 

has been especially attractive to us. Indeed, we think that the basic properties of fractional 
^ ' differential systems can be studied in a natural and nice way using algebraic integration. 

(See also \X6\, where this approach is used to study the law of the solution of a fractional 
SDE.) In the present article, we will illustrate the flexibility of the algebraic integration 
formalism by studying fractional equations with delay. More specifically, we will consider 
the following equation: 

Xt = ^0 + /o <^(-^s, Xs-n, • • • , Xs-rk)dBs + /q b{Xs, X^-n, • • • , Xs_rk)ds, t G [0, T], 
Xt = ^t, tG[-rfc,0]. 

(1) 
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Here the discrete delays satisfy < ri < . . . < < oo, the initial condition ^ is a 
function from [-rk,0] to M", the functions a : M"'^+i M"-'^, b : M"'^+i are 
regular, and S is a (i-dimensional fractional Brownian motion with Hurst parameter H > 
|. The stochastic integral in equation ([T]) is a generalized Stratonovich integral, which 
will be explained in detail in Section O Actually, in equations like ([1]), the drift term 
/q b{Xs, Xg^n, • • • , Xs-ri,)ds is usually harmless, but causes some cumbersome notations. 
Thus, for sake of simplicity, we will rather deal in the sequel with delay equations of the 
type ^ 

j Xt = ^0 + Jo cr{Xs, Xs-n, • • • , Xs-rk)dBs, t G [0, T], /^n 

\Xi = 6, tG[-rfc,0]. 
Our main result will be as follows: 

Theorem 1.1. Let ^ G C^{[-rk, 0]; R'^), a G C|(M"'^+^ R'^''^'), and let B he a d-dimensio- 
nal fractional Brownian motion with Hurst parameter H > ^. Then equation ^ admits 
a unique solution on [0, T] in the class of controlled processes (see Definition \2.5[ ) 

Stochastic delay equations driven by standard Brownian motion have been studied 
extensively (see e.g. and pj^j for an overview) and are used in many applications. 
However, delay equations driven by fractional Brownian motion have been only considered 
so far in [7j, where the one- dimensional eqyiaXioTi 

/ X, = ^0 + /o (y{Xs^r)dBs + /o h{X,)ds, t G [0, T], . . 

\X, = 6, te[-r,0], 

is studied for if > i. Observe that ([H]) is a particular case of equation 

To solve equation ([2]), one requires two main ingredients in the algebraic integration 
setting. First of all, a natural class of paths, in which the equation can be solved. Here, this 
will be the paths whose increments are controlled by the increments of B. Namely, writing 
{5z) St = Zt — Zs for the increments of an arbitrary function z, a stochastic differential 
equation driven by B should be solved in the class of paths, whose increments can be 
decomposed into 

zt~ Zs = Cs{Bt - Bs) + Psu for < s < t < T, 

with C belonging to Cl and p belonging to C^^ , for a given 7 G (|, H). (Here, denotes 
a space of /i-Holder continuous functions of i variables, see Section O) This class of 
functions will be called the class of controlled paths in the sequel. 

To solve fractional differential equations without delay, the second main tool would be 
to define the integral of a controlled path with respect to fractional Brownian motion and 
to show that the resulting process is still a controlled path. To define the integral of a 
controlled path, a double iterated integral of fractional Brownian motion, called the Levy 
area, will be required. Once the stability of the class of controlled paths under integration 
is established, the differential equation is solved by an appropriate fixed point argument. 

To solve fractional delay equations, we will have to modify this procedure. More specif- 
ically, we need a second class of paths, the class of delayed controlled paths, whose incre- 
ments can be written as 

k 

zt-zs = C^^\Bt -Bs) + Y, - Bs-n) + Psu for < s < t < T, 

i=l 



where, as above, (^^^ belongs to C7 for i = 0, . . . , k, and p belongs to for a given 
I < 7 < if. (Note that a classical controlled path is a delayed controlled path with 
(^(*) = for i = 1, . . . ,k.) For such a delayed controlled path we will then define its 
integral with respect to fractional Brownian motion. We emphasize the fact that the 
integral of a delayed controlled path is actually a classical controlled path and satisfies a 
stability property. 

To define this integral we have to introduce a delayed Levy area B^(f) of 5 for f G 
[—Tk, 0]. This process, with values in the space of matrices M*^'*^ will also be defined as an 
iterated integral: for 1 < i, j < d and < s < t < T, we set 

J s J s+v J s 

where the integral on the right hand side is a Russo-Vallois integral [21] ■ Finally, the 
fractional delay equation ([2]) will be solved by a fixed point argument. 

This article is structured as follows: Throughout the remainder of this article, we 
consider the general delay equation 

(dyt = a{yt,yt-r,,...,yt-rjdxt, tG[0,T], 

\ yt = 6, te[-rk,0], ^ > 

where x is 7-Holder continuous function with 7 > | and ^ is a 27-Holder continuous 
function. In Section 2 we recall some basic facts of the algebraic integration and in 
particular the definition of a classical controlled path, while in Section 3 we introduce the 
class of delayed controlled paths and the integral of a delayed controlled path with respect 
to its controlling rough path. Using the stability of the integral, we show the existence 
of a unique solution of equation (jlj) in the class of classical controlled paths under the 
assumption of the existence of a delayed Levy area. Finally, in Section 4 we specialize our 
results to delay equations driven by a fractional Brownian motion with Hurst parameter 
H>\. 

2. Algebraic integration and rough paths equations 

Before we consider equation (jl]), we recall the strategy introduced in [9] in order to 
solve an equation without delay, i.e., 

dyt = a{yt)dxt, t G [0, T], yo = a e M", (5) 

where a; is a R''- valued 7-Holder continuous function with 7 > |. 

2.1. Increments. Here we present the basic algebraic structures, which will allow us 
to define a pathwise integral with respect to irregular functions. For real numbers < 
a < b < T < 00, a vector space V and an integer A; > 1 we denote by Ck{[a,b];V) 
the set of functions g : [a,b]'' V such that gti - tk — whenever ti = tj+i for some 
1 < i < k — 1. Such a function will be called a (/c — l)-increment, and we will set 
C*([a, b];V) = Ufc>iCfc([a, b];V). An important operator for our purposes is given by 

fc+i 

6 : a{[a,b];V)-.Ck+i{[a,b];V), {Sgh-t,,, = Y.(-l)'-'9t....u-t,,., (6) 

i=l 



where tj means that this argument is omitted. A fundamental property of 6 is that 66 = 0, 
where 66 is considered as an operator from Ck{[a, b]; V) to Cfc+2([a, b];V). We will denote 
ZCk{[a,b];V) = Ck{[a,b];V) nKeT6 and BCk{[a,b];V) = Ck{[a,b];V) nlm6. 

Some simple examples of actions of 6 are as follows: For g & Ci{[a,b];V), h E C2{[a, b]; V) 
and / G Cs{[a, b]; V) we have 

{^9) St = 9t ~ ds, i6h)sut = hst — hsu — hut and {6f)suvt = fuvt ~ fsvt + fsut ~ fsUV 

for any s, u,v,t E [a, b]. Furthermore, it is easily checked that ZCk+i{[a, b];V) = BCk{[a, b]; 
V) for any k > 1. In particular, the following property holds: 

Lemma 2.1. Let k > 1 and h G ZCk+i{[a,b];V) . Then there exists a (non unique) 
f G Cfc([a, b]; V) such that h = 6f . 

Observe that Lemma [2.11 implies in particular that all elements h G C2{[a,b];V) with 
6h = can be written as h = 6f for some / G Ci([a, 6]; l^). Thus we have a heuristic 
interpretation of 5|c2([a,b];V)- it measures how much a given 1-increment differs from being 
an exact increment of a function, i.e., a finite difference. 

Our further discussion will mainly rely on /^-increments with k < 2. For simplicity of 
the exposition, we will assume that y = M'' in what follows, although V could be in fact 
any Banach space. We measure the size of the increments by Holder norms, which are 
defined in the following way: for / G C2([a, b]; V) let 

and 

C^{[a,b];V) = {/ G C2{[a,b];V); \\f\\, < 00}. 

Obviously, the usual Holder spaces Cf([a, 6]; V) are determined in the following way: for 
a continuous function g G Ci([a, b]; V) set 

and we will say that g G Ci{[a,b];V) iff is finite. Note that || ■ ||^ is only a semi- norm 
on Ci([a, b]]V), but we will work in general on spaces of the type 

b];V) = {g : [a, b] ^V; ga = a, Hfi-H;, < 00} , 

for a given a E V, on which is a norm. 
For h G C3([a, b]; V) we define in the same way 

\hs' t\ 

W^lkp = sup — — (7) 



:,U,tS:[a,b] 



\u - spit - u\P 



inf 1^ \\hi\\p^^^-pr^ {pi, /ii)igN with hi G C3([a, b]; V), ^hi = h,0 < pi < fi^ . 

Then || ■ ||^ is a norm on C^^la, b]; V), see [5], and we define 

C^{[a, b]; V) := {h G C^Ha, b]; V); \\h\\^ < 00} . 

Eventually, let Cl^{[a,b]]V) = U^^iCl^{[a,b]]V) and note that the same kind of norms 
can be considered on the spaces ZC^da; b]; V) , leading to the definition of the spaces 
ZC^i[a-b]-V) and ZCl+{[a,b];V). 



The crucial point in this algebraic approach to the integration of irregular paths is that 
the operator 6 can be inverted under mild smoothness assumptions. This inverse is called 
A. The proof of the following proposition may be found in [9], and in a simpler form in 

m- 

Proposition 2.2. There exists a unique linear map A : ZCl'^{[a,b]]V) — > C2^{[a,b];V) 
such that 

6A = Idzci+([a,b];v) (^^d A6 = Idci+ ^[^M-yy 

In other words, for any h G C3^([a, 6]; V") such that 6h = 0, there exists a unique g = 
A(/i) G C2'^{[a,b]]V) such that 5g = h. Furthermore, for any fi > 1, the map A is 
continuous from ZCf^{[a, h] \ V) to ([a, 6]; V) and we have 

\m\, < ^^ll/^IU, h G ZC^,{[a,h]- V). (8) 

This mapping A allows to construct a generalised Young integral: 

Corollary 2.3. For any 1-increment g G C2{[a,b]]V) such that 5g G C3^([a, 6];^) set 
Sf = {Id-M)g. Then 

n 

{6f)st = , 

list . 

1=0 

for a < s < t < b, where the limit is taken over any partition Ust = {to = s, . . . ,tn = t} 
of [s,t], whose mesh tends to zero. Thus, the 1-increment Sf is the indefinite integral of 
the 1-increment g . 

We also need some product rules for the operator 6. For this recall the following 
convention: for g G C„([a, 6]; M''"^) and h G Cm([a, K*^'^) let gh be the element of 
Cn+m-i{[a,b]-W'P) defined by 

(9) 

for ti, . . . G [a, 6]. 

Proposition 2.4. It holds: 

(i) Let g G Ci([a, 6]; M'''^) and h G Ci{[a,b],W^). Then gh G Ci(M') and 

S{gh) = 5gh + g 6h. 

(li) Let g e Ci{[a,b]-M}^'^) and h e C2{[a,b]]W^). Then gh e C2{[a,b];W) and 

S{gh) = —6g h + g 5h. 

(m) Let g G C2([a, &]; M'''^) and h G Ci{[a,b]]W^). Then gh G C2([a,6];M') and 

5{gh) = 5gh + g 6h. 

(iv) Let g G C2([a, 6]; M'''^) and h G C2{[a,b]]W^'P) . Then gh G C3([a, 6]; M''^) and 

5{gh) = —6g h + g 5h. 



2.2. Classical controlled paths (CCP). In the remainder of this article, we will use 
both the notations fdg or Jst{,f dg) for the integral of a function / with respect to a 
given function g on the interval [s,t]. Moreover, we also set ||/||oo = sup^jgi^d,! for a 

function / : M'^'' M™-'". To simplify the notation we will write Cl instead of Cj^([a, V"), 
if [a, h] and V are obvious from the context. 

Before we consider the technical details, we will make some heuristic considerations 
about the properties that the solution of equation ^ should enjoy. Set o"t = o" (yt), and 
suppose that ?/ is a solution of ([S]), which satisfies y E for a given | < k < 7. Then 
the integral form of our equation can be written as 

yt = a+ / d-udxu, t G [0,T]. (10) 
Jo 

Our approach to generalised integrals induces us to work with increments of the form 
{Sy)st = yt — ys instead of ( |TOl) . It is immediate that one can decompose the increments 
of into 

{^y)st = / cr^dxu = ^s{Sx)st + pst with pst = {^u - c^s)dxu. 

J s J s 

We thus have obtained a decomposition of y of the form 5y = a6x + p. Let us see, still 
at a heuristic level, which regularity we can expect for a and p: If a is bounded and 
continuously differentiable, we have that a is bounded and 

\^t - osl < ||o-'||oo||y|U|i - si", 

where |||/||k denotes the /t-Holder norm of y. Hence a belongs to and is bounded. 
As far as p is concerned, it should inherit both the regularities of Sa and x, provided 
that the integral J^{d'u — d's)dxu = J^{5a)sudxu is well defined. Thus, one should expect 
that p G Cl*^. In summary, we have found that a solution 6y of equation (|TOl) should be 
decomposable into 

5y = aSx + p with a E bounded and p G C'^'^. (11) 

This is precisely the structure we will demand for a possible solution of equation ([5]) 
respectively its integral form ( JTOl) : 



Definition 2.5. Let a < b < T and let z he a path in Cf([a, 6];M") with k < 7 and 
2/? + 7 > 1. We say that z is a classical controlled path based on x, if Za = a E and 
5z E ([a, b]] M") can be decomposed into 

5z = C,5x + r, i.e. {6z)st = Cs{Sx)st + Pst, s,tE[a,b], (12) 

with C G ([a, b];W''^) and p E C|''([a, b]; W). 

The space of classical controlled paths on [a,b] will be denoted by ^^.^([ct, K"), and a 
path z E QK,o([a, IR") should be considered in fact as a couple {z, (). 
The norm on QK,o([a, IR") is given by 

Af[z;Q^^a{[a,b\;R )\ = sup — — —+ sup — — + sup |G| + sup 



^,.-<^i^,^i , I s,te[a,b] \S~i^\ te[a,b] s,tG[a,b] 1^ ~ t-] 



Note that in the above definition a corresponds to a given initial condition and p can 
be understood as a regular part. Moreover, observe that a can be negative. 
Now we can sketch the strategy used in p], in order to solve equation (I5l): 

(a) Verify the stability of Qk,o([o, M") under a smooth map : R" — > R"'*^. 

(b) Define rigorously the integral J z^dxu = J{zdx) for a classical controlled path z 
and compute its decomposition ffT^ . 

(c) Solve equation ^ in the space QK,o([a, 6]; M") by a fixed point argument. 

Actually, for the second point we had to impose a priori the following hypothesis on the 
driving rough path, which is a standard assumption in the rough paths theory: 

Hypothesis 2.6. The W^-valued 'j-Holder path x admits a Levy area, i.e. a process 
= J(dxdx) e C2^{[0, T]; R'^''^), which satisfies 6x.^ = 6x ® 6x, that is 

[{6x.'^)sut] {hi) = [Sx']su[Sx^]ut, for all s,u,t e [0,T], i,j G {1, . . .,d}. 
Then, using the strategy sketched above, the following result is obtained in |9]: 

Theorem 2.7. Let x be a process satisfying Hypothesis \2.b\ and let o G C^(R"; R"'*^) he 
bounded together with its derivatives. Then we have: 

(1) Equation admits a unique solution y in ^^.^([O, T]; R") for any k < ■y such 
that 2k + 7 > 1. 

(2) The mapping (a, x, x^) ^ y is continuous from R"xC7([0, T]; R'^) xC2^([0, T]; R"'-"') 
to QK,a([0, T]; R"), in a sense which is detailed in [9, Proposition 8]. 



3. The delay equation 
In this section, we make a first step towards the solution of the delay equation 

dyt = (^iyt,yt-r^,---,yt-rjdxt, te[0,T], ^^^^ 
yt = 6, t e [-rfc,o], 

where x is a R''- valued 7-Holder continuous function with 7 > |, the function a G 
C'^(R"'^"^^; R"'"') is bounded together with its derivatives, ^ is a R"-valued 27-Holder 
continuous function, and < ri < . . . < < 00. For convenience, we set Tq = and, 
moreover, we will use the notation 

siy)t = iyt-r„...,yt-r,), tG[0,T]. (14) 

3.1. Delayed controlled paths. As in the previous section, we will first make some 
heuristic considerations about the properties of a solution: set at = a{yt,s{y)t) and 
suppose that y is a solution of ffT^ with y E for a given ^ < k < 7. Then we can write 
the integral form of our equation as 

{6y)st= / ciudxu = crs{6x)st + Pst with pst = / {an - crs)dxu. 



Thus, we have again obtained a decomposition of y of the form 6y = aSx + p. Moreover, 
it follows (still at a heuristic level) that a is bounded and satisfies 

k 



< Ik'lloo ^ \yt~r, - ys~n\ <{k+ l)||cr'||oo||l/||7|t 



1=0 



Thus, with the notation of Section [TTj we have that d belongs to CI and is bounded. The 
term p should again inherit both the regularities of and x. Thus, one should have that 
p e . In conclusion, the increment by should be decomposable into 

by = abx + p with a E Cj bounded and p G Cg*^. (15) 

This is again the structure we will ask for a possible solution to f|T3|) . However, this 
decomposition does not take into account that equation ffT^ is actually a delay equation. 
To define the integral audxu, we have to enlarge the class of functions we will work 
with, and hence we will define a delayed controlled path (hereafter DCP in short). 



Definition 3.1. Let < a < b < T and z e ([a, 6]; M") with | < k < 7. We say that z 
is a delayed controlled path based on x, if Za = a belongs to M" and if bz G ([a, 6]; K") 
can be decomposed into 

k 

{bz)st = C^^ {bx)s-ra-r. + Pst for s, t G [a, 6], (16) 
where p G C|'^([a, 6]; M") and G ([a, b]- M"'^) for i = {),... ,k. 

The space of delayed controlled paths on [a,b] will be denoted by ^'^.^([o, 1^"); and a 
path z G PK,o([a, b]] M") should be considered in fact as a {k + 2) -tuple (z, . . . , C'''^^)- 
The norm on PK,o([a, 6]; M") is given by 

Ar[^;P«,„([a,6];M")]= sup 4^^+ sup 



s,t€[a,b] |S — s,te[a,fe] P — t 



2k 



+ sup ICi^^l + sup 



iy 



St 



i=0 



t&[a,b] .^Q s,tG[a,b] \s 



Now we can sketch our strategy to solve the delay equation: 

(1) Consider the map defined on QK,a([a, b]; M") x QK,a([a — rk,b — ri]; M") by 

{T,{z,z))t = a{zt,5{z)t), tG[a,6], (17) 

where we recall that the notation s{z) has been introduced at dEl). We will show 
that To- maps QK,a([«, b]] IR") x QK,a([« — rk,b — ri]; M") smoothly onto a space of 
the form r'«,^([a,'6];R"''^). 

(2) Define rigorously the integral / Zudxu = J{zdx) for a delayed controlled path 
z G r'K,Q,([a, 6]; M""'"^), show that J{zdx) belongs to QK,a([a, 6]; M*^), and compute 
its decomposition (fT2|) . Let us point out the following important fact: T„ creates 
"delay", that is T„{z,z) G X^^.a ([«,&]; I^"'*^), while J creates "advance", that is 
J{zdx) EQ^,^{\a,b]-W). 

(3) By combining the first two points, we will solve equation (fT^ by a fixed point 
argument on the intervals [0, ri], [ri, 2ri], .... 



3.2. Action of the map T on controlled paths. The major part of this section will 
be devoted to the following two stability results: 



Proposition 3.2. Let < a < b < T. let a, a be two initial conditions in M" and let 
if e C^(M"''^+^; R') be bounded with bounded derivatives. Define on 6]; M") x 

QK,a([a -rk]h- ri]; R") by T^{z, z) = z, with 

zt = ip{zt,s{z)t), t e [a,b]. 

Then, setting a = (p{a,5{za)) = (p{a, Za-n, Za-r^_^,a), we have T^{z, z) G I^K,a([a; 6]; 
R') and it admits a decomposition of the form 

k 

{5z)^^ = Cs {5x)st + {^x)s-n,t-n + Psu s, t E [a, b], (18) 

i=l 

where (, C^*^ are the W'^^-valued paths defined by 

{zs,s{z)s),...,^ {zs,5{z)s) )Cs, se [a,b], 

OXi^O OXn,0 J 



and 

dx 

for i = 1, . . . , k. Moreover, the following estimate holds: 

Af[z;V^4[a;b];R')] (19) 
< c^,T (1 + Af^[z; Q«,„([a, 6]; R")] + M^[z; Q,,a{[a - r^^b - n]; R'^)]) , 
where the constant c^^t depends only (p and T . 
Proof. Fix s,t e [a, 6] and set 

^i'^= (^(^3,^(^)3), •••,^(^.,5(5)3) 
for i = 0, . . . , A;. It is readily checked that 

k 

= 4'^Cs{SxU + J2i^fCs-nMs-ra-n + pit + pI, 
i=l 

where 



pit - i^i'^Pst + J2^fps 



^s—ri,t—rii 
1=1 

Pst ~ '^{^t—roj ^t—m ■ ■ ■ ) ^t—rk) ¥'(^s-ro) ^s—m ■ ■ ■ ) ^s—Vk) 



s—ri,t—ri 
i=l 



(i) We first have to show that p^^p G C|'^([a, 6]; R'). For the second remainder term 
Taylor's formula yields 

\pit\ < ^ii9:'iioo(^i(fc).tr+|:i(55).-..t-..rj, 



and hence clearly, thanks to some straightforward bounds in the spaces Q, we have 

-J^. < l^'Woo (^J^'lz; Q«,„([a,6];M")] + ^Ar2[5; Q,^„([a - r„ 6 - r,]; M")] j . 

(20) 

The first term can also be bounded easily: it can be checked that 

^^fp: < llv^'lloo (^Ar[p;Cr([a,6];M")] +|jAr[p,Cr([a-r,,6-r,];M")]j (21) 

Putting together the last two inequalities, we have shown that decomposition fllSp holds, 
that is 

k 

{6z)^, = ^(°)C(fe).,, + J2 ^fCflX^x),.ra-n + Pst 

i=l 

With p,t = p\^ + pleCl%[a,b]-M.''). 

(a) Now we have to consider the "density" functions 

Clearly C, are bounded on [a, 6], because the functions t/'^*) are bounded (due to the 
boundedness of ip') and because C, C*^*-* are also bounded. In particular, it holds 

sup \Q < h'Woo sup Id, sup |Ci*)| < sup \Cs-n\ (22) 

s£[a,b] s£[a,fe] sG[a,fe] sG[a,fe] 

for i = 1, . . . ,k. Moreover, for z = 1, . . . , /c, we have 

ici?-c?i 



A: 

< W Wool Zs^ ~ Zs^\ sup ICs-rJ + llv^'lloo ^ l^si-r, - 5s2-r-J SUp |Cs-r, 
s£[a,b] s£[a,b] 



S2-ri 



<y"\\^Af[z;C^{[a,b]-R'')] sup \C-rMs2 - (23) 

sG[a,f)] 



k 



+ ll^"lloo J]Ar[5;Ci«([a-r„6-r,-];M")] sup KV.J 1^2 - 

j=l se[a,b] 

+ ||^«|U Ar[C; Ci^([a - r^, 6 - r,]; M")] - Si^- 
Similarly, we obtain 

\L-U < y"\\ooAf[^;C^{[a,b];R^)] snp \Cs\\s2-sir 

sG[a,b] 

k 

+ y'\\ooY,M[z-Ct(\a-r,,h-r,]-R^)] sup ICIka-Sir 

+ l|V^»lUAr[C;Cf([a,6];M")]|.2-sir. (24) 



Hence, the densities satisfy the conditions of Definition 13.11 

(Hi) Finally, combining the estimates fEUj) . ([21]), fl2^ and fl2^ yields the estimate ffTUj) . 
which ends the proof. □ 

We thus have proved that the map T^p is quadratically bounded in z and z. Moreover, for 
fixed z the map T<^(-, z) : I^'^) I^'')] is locally Lipschitz continuous: 

Proposition 3.3. Let the notation of Proposition \3.2\ prevail. Let < a < b < T , let 

^(1)^^(2) g Q^^^([a,6];R") and let 2 e Q^.^lja - rfc, 6 - n]; M"). Then, 

Af[Tpiz^'\z)-Tpiz^^\zy,V^,oi[a;b];R'')] (25) 
< c^,T (1 + , ^(2) ,z)y Ar[z^'^ - ; ( [a, 6] ; M")] , 

+Ar[^«; Q«,„([a,6];M")] +Ar[^(2). 5]. r")] (26) 

and the constant c^p^x depends only on ip and T . 

Proof. Denote z^^"* = for j = 1, 2. By Proposition 13.21 we have 

k 

= C?-) {6x)st + J2 (s'^ i^^)s-n,t-n + pi?, s, t e [a, b] 



i=l 



with 
where 



C?^ = Ci"^'^ = 4'''^Cs-n, s E [a, 6], 

for i = 0, . . . , /c, j = 1, 2. Furthermore, it holds = pf^^^ + p'st^\ where 

1=1 

Pst ' '~^i.^t—roi^t-rii---i^t-r^^ '^(^s— rg ' ^s-f"! ' • • • ' ■^s-^'fc ) 

Thus, we obtain for z = z^^^ — i*-^-* the decomposition 

k 
i=0 

With Ci°) = V^i^'^^d'^ - ^i°''^d'\ the paths C« are defined by ci'^ = {^"'"'^ - 4''^%-r, 
for i = 1, . . . , fc, and pst = P'st - pif • 



1;'- 'I 



In the following we will denote constants (which depend only on T and ip) by c, regard- 
less of their value. For convenience, we will also use the short notations A/'[z], A/'[z*-^-*], 
Mlz^"^^] and J\f[z^^^ — z^"^^] instead of the corresponding quantities in ( |25l) - (l26l) . 

(i) We first control the supremum of the density functions C i = 0, . . . ,k. For i = 0, 
we can write 

= ^rnc« - ) + (^r^ - ^r^)^) 

and thus it follows 

<c(l+Ar[^(2)]) _^[^(i)_^(2)] (27) 

Similarly, we get 

\C^^\<cAr[z]Ar[z^^^ -z^% (28) 

(ii) Now, consider the increments of the density functions. Here, the key is to expand the 
expression ipi^'^^ — ipi^''^^ for i = 0, . . . , k. For this define 

^«-zf) + zf, rG[0,l], se[a,b]. 
'A'\<~^)s) = ^(«.(l),5(z).)-^(«.(0),5(z).) 

Hence it follows 

Note that 6**^''*-' is clearly bounded and, under the assumption Lp G Cf , it moreover satisfies: 

_ < c (Ar[^(i)] + Ar[^(2)] 1^ _ ^|._ (3Q) 

For i = we can now write 

go) _ 4(0) ^ ^^(0,1) _ ^(0,1) j (^1) _ ^2)) ^ ^(0,) ^(gl) _ g2)) _ (^l) _ 4,))^J 

+ (^f'^^ - icP - CP) + cf ^ ((^f - ^f'^^) - - 

It follows 

_ (?(o)| < c (Ar[2(i)] +Af[z]) \t - A/'[^(i) - ^(2)] + cAr[z(i) - \t - 

+Ar[.(^)] (#^)-#^V(^r)-^r^) . (31) 





Us[r) = r[z 


We have 




d(p 








where 





s 



Using (I22D and 1^ we obtain 

< c (1 + Ar[zW] + Ar[2(2)] + Ar[z« - ^(2)] |t _ 

Combining flHTj) and (13^ yields 

|go) _ ^(o)| < ^ (1 +^[^(1)] +_^[^(2)] + - ^(2)] 1^ _ 

By similar calculations we also have 

ICP - Ci'^l < c (1 + Ar[^«] + Ar[z(2)] + Ar[z])' Ar[^« - |t - 
for i = 1, . . . , /c. 

(iii) Now, we have to control the remainder term p. For this we decompose p as 

A _ .(1) , .(2) 

Pst — Pst + Pst ! 

where 

1=1 



i=l 

We consider first p*^^); for this term, some straightforward calculations yield 

Ipii^l < c{l+Af[z^^^]+Af[z])Af[z^'^ - \t - s\^\ 
Now consider p^'^\ The mean value theorem yields 

pi? = - V^i^'^O (5^^^^).. - m'''' - i^r') {Sz^'')st 

k 



It'' ' % 
i=\ 



(^(0,1) _ ^(0,1)) (5(^(1) _ ^(2))) 

+ ((^(0,1) _ ^(0.2)) _ (^(0,1) _ ^(0.2))) (5,(2) 



St 

M) „;,(o,2)^^ ^x.(2)\ 

St 



i=l 



^Qi + Q2 + Q 



3) 



with 

We shall now bound Qi, Q2 and Q3 separately: it is readily checked that 

<c(i+Ar[2W]+A/-[2(2)]+Ar[5]) \t-s\\ 

and thus we obtain 

In order to estimate Q2 and Qs, recall that by ( |29l) in part (^m^ we have 

^(.1) _ ^(.2) ^ _ ^(2))^ _ _ _ ^ _ ^(2))) ^ ^37) 

where 

M,(r')=;2« + r'(zp)-z«). 
Similarly, we also obtain that 

^(.1) _ ^(.2) ^ _ ,(2)), . . . , _ ,(2))) (38) 

with 

Now, using (lH7j) and (IHHl) we can write 

for any i = 0, . . . , k. Since moreover 
Us{r,r') - {us{r'),5{z)s) 

= r ^(zj ^ — 2^ ^) + r (^j^ ^ — ^ — (zj ^ — 2^ ^)), — . . . , Zf^Vk " ^s-r^j , 

another Taylor expansion yields 

|^(M) _^(/,*))| <c (Ar[z«]+Ar[z(2)]+Ar[5]) \t-s\\ 

Hence, we obtain 

< c (Ar[2(^)] + Af[z^^^] + Ar[5]) Ar[^(i) - |t _ (39) 

from which suitable bounds for Q2 and Q3 are easily deduced. Thus it follows by ( l36l) 
and ([39]) that 



Combining this estimate with flHoj) we finally have 

\pst\ < c (1 + Af[/^^] + Ar[z(2)] + _^f[z]Y Af[z^^^ - |t _ s|2-. (40) 

(iv) The assertion follows now from ^ and 

□ 

3.3. Integration of delayed controlled paths (DCP). The aim of this section is to 
define the integral J'{m*dx), where m is a delayed controlled path m G I^*^)- 
Here we denote by A* the transposition of a vector or matrix A and by Ai ■ A2 the inner 
product of two vectors or two matrices Ai and A2. We will also write Q^.a (resp. V^^a) 
instead of QK,a([«, V) (resp. V^^a{[0',b]', V)) if there is no risk of confusion about [a,b] 
and V. 

Note that if the increments of m can be expressed like in (fTGll . m* admits the decom- 
position 

k 

(Sm*),, = + Pit, (41) 

i=0 

where p* E C|'*([a, 6]; M"'^''^) and the densities i = 0, . . . , /c satisfy the conditions of 
Definition 13.11 

To illustrate the structure of the integral of a DCP, we first assume that the paths 
C^'^ and p are smooth, and we express J'{m*dx) in terms of the operators 6 and A. In 
this case, J{m*dx) is well defined, and we have 

mldxu = m*{xt ~ + ^ ("^« " ^*s)d^u 

for a < s < t < 6, or in other words 

J{m* dx) =m* 5x + J {6m* dx). (42) 
Now consider the term J'{Sm* dx): Using the decomposition fl41l) we obtain 



k 



J{6m* dx) = / f E (^^):-n,«-nCi^^* + Plu] dxu = A^t + Jstip* dx) (43) 



with 



Since, for the moment, we are dealing with smooth paths, the density can be taken 
out of the integral above, and we have 

k 

with the d X d matrix ^^t{v) defined by 



= { I { I dx^j dx^u ,•••,/ I / dx,„ \ dx^^ ] , < s < t < T 

s-\-v / J s \J s-\-v 



for V G {— T/c, . . . , — ro}. Indeed, we can write 

J s J s 

J s 

Inserting the expression of Ast into (H2l) and ( H3l) we obtain 

J7.i(m* dx) = m:{Sx),t + Ci'^ " ^li-r^) + Jst{p* dx) (44) 
for a < s <t <h. 

Let us now consider the Levy area term x^^(— rj). If a; is a smooth path, it is readily 
checked that 

[5^^{-ri)]sut = ^li-n) - x^„(-r,) - x^i(-ri) = {6x)s-r,,u-n ® {Sx)ut, 

for any i = 0, . . . , k. This decomposition of 5x^(— rj) into a product of increments is the 
fundamental algebraic property we will use to extend the above integral to non-smooth 
paths. Hence, we will need the following assumption: 

Hypothesis 3.4. The path x is a W^-valued •y-Holder continuous function with 7 > | 
and admits a delayed Levy area, i.e., for all v G {— r^, . . . , — ro}, there exists a path 
^^{v) G C2^([0,T];M'^''^), which satisfies 

5^^{v)=5x'' ®5x, (45) 

that is 

[{5y?{v))sut] = [Sx']s+v,u+v[Sx^]ut for all s,u,t E [0,T], i,j G {l,...,d}. 
In the above formulae, we have set x"" for the shifted path x^ = Xg+v 

To finish the analysis of the smooth case it remains to find a suitable expression for 
J^{p* dx). For this, we write (H^ as 

k 

Jstip* dx) = J,t{m* dx) - m:{6x),t - J2 ^ " ^'*(-^^) (^6) 
and we apply 6 to both sides of the above equation. For smooth paths m and x we have 



5{J{rn* dx)) = 0, 6{m* 6x) = —6m* 6x, 



by Proposition 12.41 Hence, applying these relations to the right hand side of fHUj) . using 
the decomposition f HTl) and again Proposition 12.41 we obtain 

[5{J{p* dx))]sut 



fc 

1=0 

k k 



i=0 i=0 



i=0 

In summary, we have derived the representation 

k 

S[J{p* dx)] =p*6x + J2 '^C^'^ ■ x'(-rO, 

i=0 

for two regular paths m and x. 

If m, X, i = 0, . . . , /c and are smooth enough, we have 5[J'{p* dx)] E ZC]^ and 
thus belongs to the domain of A due to Proposition 12. 2[ (Recall that = 0.) Hence, it 
follows 



J{p* dx) = a(^p*6x + J2 ^(^^^ ■ x^(-^i) j 



and inserting this identity into (jS]), we end up with 

k / k \ 

J{m* dx) = m*6x + ^ C^'^ ■ x^(-ri) + A p*6x + ^ 6('^'^ ■ ^^(-r,) J . (47) 

The expression above can be generalised to the non- smooth case, since J{m* dx) has 
been expressed only in terms of increments of m and x. Consequently, we will use ( H71) as 
the definition for our extended integral. 

Proposition 3.5. For fixed ^ < k, < 'y, let x be a path satisfying Hypothesis Fur- 
thermore, let m G 'DK,a([a, 6]; M*^) such that the increments of m are given by fTB) . Define 
z by Za = a with a G M and 

k f ^ \ 

i=0 \ i=0 j 

for a < s < t < b. Finally, set 

J{m* dx) = 6z. (49) 

Then: 

(1) J{m* dx) coincides with the usual Riemann integral, whenever m andx are smooth 
functions. 



(2) z is well-defined as an element of QK,o([a, b]; M) with decomposition 6z = ■m*6x + p, 
where p G C|''([a, b]; M) is given by 



P 



J2 C^'^ ■ x'(-r.) + A p*5x + J2 SC^'' ■ x'(-rO 



1=0 



(3) The semi-norm of z can he estimated as 

M\z- Q,,„([a, b\- M)] < ||m|U + c.„t(6 - a)^-'^Ar[m; P«,^([a, b\X)\ 
where 



— r 



i)\\2j 



i=0 



with the constant CK,7,<^,r depending only on Kjjjip and T. Moreover, 



(4) It holds 

J'stim* dx) = lim > 



i=0 



(50) 



(51) 



(^^).,...+Eci/^-xl,,,(-r,) (52) 

for any a < s < t < b, where the limit is taken over all partitions list = {s = 
tQ, . . . ,tN = t} of [s, t], as the mesh of the partition goes to zero. 



Proof. (1) The first of our claims is a direct consequence of tlie derivation of equation 

dSD- 

(2) Set Cx = + X]i=o I|x^(^'''i)ll27- Now we show that equation (HSjl defines a classical 
controlled path. Actually, the term m* Sx is trivially of the desired form for an element 
of Qk,o- So consider the term h^^^ = XliLoCi*^ ' ^^ti^^i) for a < s < t < b. We have 

\h^st^\ < E IIC«l|oo|x?,(-r.)| < IIC^Iloo) c.|t - s\'- 

i=0 \ i=0 J 



Thus 



The term 



\\h'^'m2.<c,{b-ay'^'''''W[m-V,Ahb];R% 



satisfies Sh^"^^ = 0. Indeed, we can write 

k 
i=0 



by Proposition 12.41 and because 66 = 0. Applying fH31) to the right hand side of the above 
equation it follows that 

k k 

6h^^^ = 6p*6x + J2 ^C^'^ ■ {Sx-''^ ® 6x) = 6p*6x + ^ 6x-''^*6C^'^*6x 

k 

= {6p* + Y,^X~"'*^C^"^*)^X■ 



^=0 



i=0 



1=0 



However, due to Proposition 12.41 it holds 



i=0 1=0 

Since the increments of m are given by ( JT6l) we finally obtain that 

5/^(2) = 6{6m*)6x = 0. 
Moreover, recalling the notation ([7]), it holds 

\\p*6x\\2fi,K < Cx{h - ay'''\\p\\2n 

and 

k k 

II 5^5C^^^ ■ x^(-r,)|U,2« < c.(6 - af^^'~^^ ^ ||C«|U. 



i=0 



i=0 



Since 7 > K > I and 6h^'^^ = 0, we have h^'^^ E Dom(A) and 

||/^('^||3«<c.(l + T^-'^)(6-ar-WKP«,^([a,6];M'^)]. 
By Proposition 12.21 it follows 



|A(/^^^^)I|3.<^I|/^(^)||3. 



and we finally obtain 



o3k 1 

II/.W - A(/.(2))||2, < c.___(i + T^-)(6_a)^-WKP.„ 



(53) 



Thus we have proved that p G Cf^da, b]; M) and hence that z G QK,a{[C', b];M). 

(3) Because of {6z)st = m*{6x)st + Pst and m G "D^.tilict, the estimates (1501) and ( 151 

now follow from fl53l). 



(4) By Proposition 12.41 (^nj and the decomposition ffTUl) we have that 

k 

6{m*6x)sut = -{6m*)su{6x)ut = -pl^{6x)ut - ^{Sx)l_,^.^^^_^^Cl'^*{6x)ut- 



i=0 



Thus, applying again Proposition 12.41 (ii), and recalling Hypothesis 13.41 for the Levy area, 
we obtain that 



6lm*6x + J2 C^'^ ■ 



i=0 



i=0 



Hence, equation fHHj) can also be written as 

Jim* dx) = [Id - A6] (m* 6x + J2 ' 



i=0 



and a direct application of Corollary 12.31 yields (l52l) , which ends our proof. 

□ 

Recall that the notation A* stands for the transpose of a matrix A. Moreover, in the 
sequel, we will denote by Cnorm a constant, which depends only on the chosen norm of 
M"''^. Then, for a matrix- valued delayed controlled path m e I^re,a([a, b]; M"''^), the integral 
J{mdx) will be defined by 



J{mdx) = {J{m^^>dx),...,j{m^''>dx))' 



1 



where m*^*^ G T>i^^s,{\a^ 6]; W^) for i = 1, . . . , n and we have set m = {m'^^\ . . . , m*^"))*. Then 
we have by (Hlj) that 

Ar[:r(mdx);Q,,„([a,6];M")] (54) 

For two paths m^^\m^'^'> G V i^^s,{\a^h\]W^'''') we obtain the following estimate for the 
difference of z^^^ = J'(m^^^ dx) and z^'^^ = J'{m^'^^ dx): As above, we have clearly 

Ar[z« - ^(2). Q,^o([a,6];R")] < c™||m« - m^^)]]^ 

+ c„ c„i(6 - a)^-^Ar [m(i) - m(2); I),_o([a, b]; R^'^)] . 

However, since m^^ = m^p it follows 

Af[z^'^ - z('); Q«,o([a, 6]; M")] < 2c™ c^b - af-'^N [m^^^ - m^^); P,_o([a, 6]; M"-'^)] . 

(55) 



4. Solution to the delay equation 

With the preparations of the last section, we can now solve the equation 

dyt = (T{yt, yt-n, • • • , yt-rj dxt, t g [0, T], 
yt = ^t, tG[-r,0], 

in the class of classical controlled paths. For this, it will be crucial to use mappings of 
the type 

r : Q,,,([a, b]; R") x Q,,^([a -n,b- n]; R'^) ^ Q,,,([a, b]; R") 

for < a < b < T, which are defined by {z, z) i— > z, where zq = a and Sz given by 
6z = J'{Tcr{z, z)), with Tcr defined in Proposition 13. 2[ /^From now on, we will use the 
convention that Zt = Zt = Zt = for t G [— r, 0]. Note that this convention is consistent 
with the definition of a classical controlled path, see Definition 12.51 since ^ is 27-Holder 
continuous, it can be considered as a part of the remainder term p. 



The first part of the current section will be devoted to the study of the map T. By 
we have that 

Ar[J(T.(z,5));Q,,„([a,6];M")] 

< c„ (Iklloo + cUb - ay-^Af[T^{z, z); V,^&{[a, 6]; M"''^)]) . 

Since 

T^{z,z) = {T^w{z, z), . . . ,T^(n){z, z))* , 
where a^^) G C3(M«'^; RL'^) for i = 1, . . . , n and a = {a^^\ . . . , a^"))*, it follows by (USD 
that 

Af%{z, z);I).,A([a; 6];M")] 

< c.,T (1 + Ar2[^; Q,,„([a, 6]; M")] + Ar2[5; ^.^^([a - r^, 6 - n]; M")]) . 
Combining these two estimates we obtain 

Ar[r(^,z);Q,,,([a;6];M")] (57) 
< Cgron,tk (1 + Af'[z; Q.A[a -rk,b- n];R'% {l + {b - ay^Af^lz; Q,,„([a, 6]; M")]) , 

where the constant Cgrowth depends only on Cint, Cnorm, cr, n, 7 and T. Thus the semi-norm 
of the mapping F is quadratically bounded in terms of the semi-norm of z and z. 

Now let z^^\z^^'^ e QnAh b]] K") and z e 2^,5 ([a -rk,b- ri];R'^). Then, by (ESI) we 
have 

Af[T{z('\S) - T{z('\~zy, Q.,o([a, 6]; M")] (58) 

< 2c™c,„i(6-a)^-W[T.(z«,z) -T.(z(2\z);P,,o([a,fe];M"'^)]. 

Applying Proposition l3.3l i.e. inequality (^3]) . to the right hand side of the above equation 
we obtain that 

AA[r(z«, ~Z) - r(^(2)^ ^n)] ^59) 

< cup{l + C{z^'\ z^'\z)Y^[z^''> - Q^^o([a, 6]; M"'^)] (6 - a)^-^ 

with a constant C;jp depending only on Cj„f, c„orm, c", ^, 7 and T, and moreover 

Ciz^'\ z^^\z) = M[z- Q.,5([a - r,, 6 - n]; M")] 

+ Ar[z«; Q.,„([a,6];M")] +Ar[z(2); Q,_„([a, 6]; M")]. 

Thus, for fixed z the mappings r(-, 5) are locally Lipschitz continuous with respect to the 
semi-norm ATf-; Qk,o(['2, b]] R")]. 

We also need the following Lemma, which can be shown by straightforward calculations: 

Lemma 4.1. Let c, a > 0, r G [0,T] and define the set 

= {m G : c(l + rV) < u}. 
Set also T* = (8c2)-V". Then we have ^'f ^ and sup{M ; u G ^^f } < (4 + 2V2)c. 

Now we can state and prove our main result: 

Theorem 4.2. Let x be a path satisfying Hypothesis \3.4\ let ^ G Cf'^{[—r,0];W^) and let 
a G CKM^'^+^M"'"^). Then we have: 



(1) Equation 13^) admits a unique solution y in QK,go([0, T]; M") for any | < k < 7 
and any T > 0. 

(2) Let F : C^^\[-r,0];W) X CJ{[0,T];R'^) X {Cl\[0,T];R'^''^)y^^ ^ C^{[0,T];W) be 
the mapping defined by 

F (e, X, x2(0), x2(-ri), . . . , x2(-rfc)) = y, 

where y is the unique solution of equation fSU) . This mapping is locally Lipschitz 
continuous in the following sense: Let x be another driving rough path with cor- 
responding delayed Levy area x.^{—v), v G {—rk, . . . , — tq}, and ^ another initial 
condition. Moreover denote by y the unique solution of the corresponding delay 
equation. Then, for every N > 0, there exists a constant Kjq > such that 



\y ?/||re,oo 
^ I ||x — 3;||7,oo 



k \ 

+ J2^f[^'{-r■^ - i'(-n); Cl^iiO, T]; R'^)] + U- ^2^,00 ] 
1=0 / 



holds for all tuples x, x^, x^(— ri), . . . , x^(— r^)), x, x^, x^(— ri), . . . , x^(— r^)) 
with 

Ar[x^(-r,); ([0, T]; R'^)] + f^^iSt^-r,); ([0, T]; R'^)] 

1=0 i=0 

+ \\x\\^,oo + ||5|l7,oo + 11^1127,00 + ||f||27,oo < N, 

where ||/||^,oo = ll/lloo + l^/U denotes the usual Holder norm of a path f . 



Proof. The proof of Theorem 14.21 is obtained by means of a fixed point argument, based 
on the map T defined above. 

1) Existence and uniqueness. Without loss of generahty assume that T = Nri. We will 
construct the solution of equation fl3Ul) by induction over the intervals [0, ri], [0, 2ri], . . ., 
[0, Nri], where we recall that ri is the smallest delay in fl3Bl) . 

(i) We will first show that equation (!56l) has a solution on the interval [0,ri]. For this 
define 

n = (8c?)-i/(^-'^)Ari, 

where 

ci = Cgro^th (1 + ^f''[^■, C'2\[-rk, 0]; R")]) . 
Moreover, choose ti G [0, fi] and Ni eN such that NiTi = ri, and define 

= [{i - '^)ri,in], i = l,...,Ni. 

Finally, consider the following mapping: Let Fi^i : QK,$o(-^i,ij Qk.^oI-^i.Ij given 

hj z = Fi 1(2;), where 

{Sz)st = Jst{Ta{z,C) dx) 

for < s < t < ri . 



Clearly, if z^^'^^ is a fixed point of the map i), then 2;'^^'^^ solves equation (1561) on the 
interval Ji^i. We shall thus prove that such a fixed point exists. First, due to fl57|) we have 
the estimate 

^f[rlAz)■ Q.,^„(/i,i;M")] < ci {l + n^-''Af'[z; M")]) . (60) 

Thanks to our choice of Ti and Lemma [4.11 we can now choose Mi G A'^i''^~^ accordingly 
and obtain that the ball 

Bm, = {ze Af[z; M")] < Mi} (61) 

is left invariant under Fi^i. Now, by changing ri to a smaller value (and then A'^i accord- 
ingly) if necessary, observe that Fi^i also is a contraction on -Bmi, see ( l59l) . Thus, the 
Banach theorem implies that the mapping Fi^i has a fixed point, which leads to a unique 
solution z^^'^^ of equation ( |56l) on the interval Ji^i. 

If Ti = ri, the first step of the proof is finished. Otherwise, define the mapping F2,i : 
Q^^^(M) (/2,i; K") Q,,4M) (/2,i; K") by i = F2,i(2;) where 

(5^)st = Jst{Taiz,^) dx) 

for Ti < s < t < 2ri . Since ri < ri, it still holds 

Ar[F2,i(z);Q^^4ia)(/2,i;M")] <ci(l + ri^-W2[z;Q^_4M)(/2,i;K")]) (62) 

and we obtain by the same fixed point argument as above, the existence of a unique 
solution 2;^^'^) of equation on the interval /2,i. 

Repeating this step as often as necessary, which is possible since the estimates on the 
norms of the mappings Fj 1, j = 1, . . . , A^i are of the same type as (137|) . i.e. the constant Ci 
does not change, we obtain that z = Xlj^i ^^'''^^ ^ij 1 the unique solution to the equation 
( l56l) on the interval [0,ri]. 

Now, it remains to verify that z given as above is in fact a CCP. First note that by 
construction z is continuous on [0, ri] and moreover that z is a CCP on the subintervals 
Ij^i with decomposition 

(52),i = CP(5x),, + pif\ s,tEl,,u 
for s <t. Clearly, we have 



{Sz)st = J2^Sz'^^'^^),,s/t„tAt,+,, s,t e [0,ri], 



J=Js 



for s <t, where tj = (j — l)ri and js,jt £ {1, • • • , ^ — 1} are such that 
Setting 



tj, < S < tj^+i < . . . < tj^ < t < tji+i. 



JVi 



and 



Jt Jt 

J=Js J=Js 



we obtain 

{Sz)st = Cs{Sx)st + Pst, s,t e [0, n] 

for s < t. 

Now, it follows easily by the subadditivity of the Holder norms that 

Ni 



and 



sup — < > sup — — 

MG[0,ri] \s-t\^ ~[s,t&I,,, \s-t\^ 



sup IGI = sup sup IQ 'I, sup -rrz<y^ sup . 

tG[0,Ti] i=l,...,Ar-ltG/j,i s,tg[0,ri] |S-t|'" ~[ s,t(^I,,i \S - t\'^ 

Furthermore, we obtain 

Thus, we have in fact that z G Q^.^odO? ^")- 

('zij Let / = 1, . . . , — 1 assume that z G Qk,^o([0) ^'^i]) ^"') is the solution of the delay 
equation ( l56l) on the interval [0,/ri]. Now we will construct the solution on the interval 
[/ri,(/ + l)ri]. Set 

Q+i = cgrowth (1 + Af^[z; Qn,zir,-r,&n - rfc, In]; M"]) 

and define 

Furthermore, choose r;+i G [0,f;+i] and A^^^+i G N such that Ni^iTi^i = ri, and define 

Ji,i+i = [In + {i- 1)t/+i, /ri + iTi+i], i = 1,..., Ni+i. 

Consider the mapping Fi^z+i : Q^.i,,^ M") Q^.i,,.^ K") hy z = Fi,/+i(z) 

where 

(5^)st = Jst{Ta{z,z) dx) 

for /ri < s < t < /ri + r^+i . Again ^(1''+^) is a fixed point of the map Fi if and only if 
solves equation (1561) on the interval However, by fl57p we have the estimate 

Ar[Fi,;+i(z); Q«,.-,,^(/i,/+i;M")] < Q+i {l + ri+{'-^M^[z; Q.,^,^ M")]). 

Now we can apply the same fixed point argument as in step (i), which leads to a unique 
solution of ( !56l) on the interval h,i+i- 

If r^+i 7^ ri, define for the next interval l2,i+i the mapping 



F2J+1 : Q (i,i+i)(/2,«+i;K") ^ Q„ ^(i.i+i)(/2,, 



hj z = F2,i+i(2;), where {6z)st = Jst{Ta{z, z) dx) for In + ri+i < s < t < In + '^n+i ■ 
Since hi + r^+i < (/ + l)ri, we still have the estimate 

Ar[F2,mW; Q^,(i.m)(/2,m;M")] < q+i (1 + r;'-,''M^[z; Q^^(M+i)(/2,m;M")]). 

Now the existence of a unique solution of (l56l) on the interval /2,i+i follows again 

by the same fixed point argument. 



Proceeding completely analogous to step (i)we obtain the existence of a unique path z G 
Qk,5,^^ ([^'"i, I^"), which solves the delay equation ( l56l) on the interval [/ri, (/+l)ri] 

for a given "initial path" z G Qk,^(,([0, /ri]; M"). Patching these two paths together, we 
obtain (using the same arguments as at the end of step (i)) a path z G Qk,5o([0' + 
l)ri]; R"-), which solves equation fl56l) on the interval [0, (/ + 

Thus we have shown that there exists a unique path z G Q^^(^^{[Q,T\,W^), which is a 
solution of the equation fl^ . Moreover, by the above construction we obtain the following 
bound on the norm of this path: 

Ar[^;Q.,5„([0,T];M")] (63) 

< / (aA[x;C7([0,T];R")] + ^Ar[x^(-r,);C2'^([0,T];M")] + A^i^; C^([0, T]; M")] 

where / : [0, oo) — > (0, oo) is a continuous non-decreasing function, which depends only 
on K, 7, n, d, cr, T and ri, . . . , r^. 

Continuity of the ltd map. Let y = F {C,, x, x^(0), x^(— ri), . . . , x^(— r^)) . Since y solves 
equation we have (5|/)st = >Jst{o'{ys,s{y)) dxs). It follows by the Propositions 13.21 and 
13.51 that 

k f ^ \ 

{5y)st = ms{6x)st + Yl Ci*^ " ^st(-r^) + A.t [pSx + J2 ^C^'^ ■ x'(-r,) (64) 

j=0 \ i=0 / 

for < s < t < T, with 

m, = a(?/„s(?/),), Cf^ = ^i'^m^-r^, ^jjf = { -^{y,,s{y)s), • • • , -^{,ys.5{y)s 



(65) 

for i = 0, . . . ,k. Moreover, note that the remainder term p of the decomposition of y 
satisfies the relation 

k ^ ^ \ 

Pst = J2 Ci^^ ■ + [pSx + Y, ^C^^^ ■ x^(-r.) . (66) 

i=0 \ i=0 ) 

Now consider fl3U|) with a different initial path ^, driving rough path x and corresponding 
delayed Levy area x^(f ), for v G {— Tfc, . . . , — ro}. If the assumptions of the theorem are 
satisfied, then also the equation 

dyt = a{yt,yt-n,---,yt-r,,)dxt, t E [0,T], 

yt =L te [-r, 0] 

admits a unique solution y = F(,^, 5(0), x^, x^(— ri), . . . ,x^(— r^)). Clearly we also have 
in this case 

k f \ 

= m,{Sx),, + J2 Ci'^ ■ (x'(-r.))si + Ast [pSS: + J2 ^C^'^ " ^'(-^^) (67) 

i=0 \ i=0 / 

for < s <t <T, with m, and defined according to ( 165|) and ( l66l) . 
(^zj We first analyse the difference between p and p. Here we have 

Pst-Pst = e'^st +Ast{e^^^), (68) 



with 



i=Q 
2/ 



e(2) = + .x2^_^,-j_ 



1=0 



1=0 



Now set 



C{y) = ||x||oo + Ikll^ + J]||x2(-r,)||2^ + A^[y;Q«,„([0,T];M-)] + ||e|U + llClh^, 



1=0 



define C{y) accordingly for y, and let R be the quantity 



R ^ \\X- X\\^ + \\X - X\\^ + J2 - X^(-^0I|27 + U - CIloo + lie - CII27- 

1=0 

In the following we will denote constants, which depend only on k, 7, n, d, a and T, by c 
regardless of their value. 

Fix an interval [a,b] C [0,T]. By straightforward calculations we have 



|e«| < c (1 + - s\'^ R + cC{y)\t - s\'^ sup - 

Te[{s-rk)+,t] 

for s, t e [a, 6]. Now, consider the term e^^^. We have 

k k 

i=0 i=0 
= (P - p)su{Sx)ut + Psu{S{x - X))ut 

1=0 1=0 

for s,u,t e [a, 6]. Clearly, it holds 

|(p - P).„(5a;)„,| < Ciy)\t - u\''\s - u\^^M[p - ~p; Ci^{[a, b];W% 
\Psu{S{x - x)U < \t-u\^\s-u\^^M[p;C^^{[a,b];R'')]R 



(69) 



and 



5^(<5(C«-C«))..-x^,(-n) 



i=Q 



< 



cC(y)|t-i.r^5;]|(5(C«-C«)). 



1=0 



j=0 



Furthermore, we also have, for any i = 0, . . . , A; that 



Recall that the Holder norm of a path / is defined by 

U,oo,[s,t] = sup \fr\+ sup 



l/n 




k2 


- Ti ^ 



Set also C = c(l + C{y) + C{y)), where c is again an arbitrary constant depending only 
on K, 7, n, d, a and T. Using these notations and combining the previous estimates, we 
end up with: 



.(2) 



< C\t - u\^\s - ul^'^R + Clt - u\^^\s - mI"^ 



\K,co,[a-ri,b-ri\ 



1=0 



+ C\t - u\'\s - uf'^MXp - p; ([a, 6]; M")]. (70) 
Hence e'-^-* belongs to Dom(A) and we obtain by Proposition 12.21 that 

k 

||A(e(2))||3. <CR^CY,\\y- y|U,oo,[a-n,6-n] + C Ar[p - p; Cl^{\a, h\-W)\. (71) 

i=0 

Inserting the estimates for e*^^^ and A(e*^^^), i.e. (lU^ and (I7T|) . into the definition flU^ of 
p — p gives finally 



Ar[p - p; ([a, 6]; M")] < C|6 - aX<-^ ||y - y|U,oo,[a-n,6-n] 

i=0 

+C|6 - ap-" i? + C|6 - ap-" Ar[p - p; C2^'^([a, 6]; R")], 
and due to the subadditivity of the Holder norms, we get 

Ar[p-p; ClWa, h\- R")] < C\h-a\^~- M\p-~p; Cl^(\a, h\-W)\^C\h-a\'-'^ ||y-y|U,oo,M] ^ 

+ C\h - a\<-^ \\y - y\U^Xa~r,A R^C\h- a\<-^ R. (72) 

(a) Now consider the difference between y and y. Completely analogous to step (i) we 
also obtain that 

Af[y-y; C^i[a, b]; R")] < C\b~ar^ Af[p- p; Cl^i[a, b]; R-)] + C\b-ar^ \\y-y\Uoo,[aM R 

+ C\b - a\''~^ \\y - y\\.,oo,[a-r„a] R + C\b- a\^'-^ R. 



Moreover, since 



sup \yr-yr\ < \ya-ya\ + {b-arAf[y-y;C^{[a,b];R^)], 

T&[a,b] 



we also have 



\\y - vUaM < C\b - a\^-^Af[p - p; Cl^{[a, 6]; M")] + C\b - a\^-^ \\y - y\\.,oo,laM R 

+ C\b - a\^-^ \\y - y\\.,oo,[a-r„a] R+\ya- Val + C\b - a\^-^ R. (73) 



(in) Now set 



^(aM=J^[p-p-A\[a,b]-R-)] + 



I K,oo,[a,b] • 



By combining fl72|) and flTHj) we finally have that 

A(a, b) < C(l + - ap-'' A(a, 6) + C(l + - ap"^ A((a - rfc) + , a) 

+ \ya-ya\ + C{l + R)\b-a\^-^R. (74) 

/ , \l/(7-«) 

Now choose a = and 6i = ( 2c(i+r) ) ■ ■'■^ ^'^^^ case, we obtain from ( 1741) that 

A(0,6i)<^A(0,6i) + |eo-eo| + ^i?, 

which yields 

A(0,6i) </? + 2|eo-eo| <3i?. (75) 
For the next interval [hi, 2bi], we obtain in turn that 

A(6i, 2b,) < iA(6i, 26i) + ^A(0, 6i) + \y,, - J + ^i?, 

and hence 

A(6i,26i) < A(0,6i)+2||/,, -yfcj+i?<10/?, 

by dZS]). 



Repeating this step [T/6iJ-times we obtain that there exists a continuous non-decrea- 
sing function g : (0, oo) (0, oo) such that 

A{tb,,{t + l)b,)<giT/b,)R 

for all i = 0, . . . , [T/bi\. Using the subadditivity of the Holder norms, we obtain the 
estimate 

A{0,T)<{l + T/bi)g{T/b,)R. (76) 
Now recall that C = c{l + C{y) + C{y)) and note that R < c{C {y) + C {y)) . Thus we have 
T/6i = T (2C(1 + i?))'/(^-") < c{C{y) + C{y))fl^^^^\ 

where 

k 

C{y) = llxlloo + + J]||x2(-r,)||2, + Ar[y;Q.,,([0,T];M")] + ||e||oo + IKib,, 

and C{y) is defined accordingly. However, by (lU^ it follows that 

C(y) + C(y)<D + /(D) + D + /(D), 

where 

^ = ||a:||oo + + ||x2(-r,) ||2^ + + H^lh^, 

and -D is again defined accordingly. Thus, we obtain now from flTHl) that there exists 
a continuous function g : [0, oo) ^ [0, oo), which depends only on K,'y,a,n,d,T and 
ri, . . . , Tfc, such that 

A(0,T) < g{D + D)R. 

Hence, the assertion follows. 

□ 



5. Application to the fractional Brownian motion 



All the previous constructions rely on the specific assumptions we have made on the 
path X. In this section, we will show how our results can be applied to the fractional 
Brownian motion. 

5.1. Definition. We consider in this section a rf-dimensional fBm with Hurst parameter 
H defined on the real line, that is a centered Gaussian process 

B = {Bt={Bl...,Bt);teW.], 

where B^,...,B'^ are d independent one-dimensional fBm, i.e., each 5* is a centered 
Gaussian process with continuous sample paths and covariance function 

RH{t,s) = l{\sr+\tr-\t-sr) (7?) 

for i = 1, . . . ,d. The fBm verifies the following two important properties: 

(scaling) For any c> 0, B^"^ = c^B./^ is a fBm, (78) 

(stationarity) For any h eW, B.^h — Bh is a fBm. (79) 
Notice that, for Malliavin calculus purposes, we shall assume in the sequel that B is 
defined a complete probability space (fi, J-', P), and that JF = a{Bs; s G M). Observe also 
that we work with a fBm indexed by M for sake of simplicity, since this allows some more 
elegant calculations for the definition of the delayed Levy area. 

5.2. Malliavin calculus with respect to fBm. Let us give a few facts about the 
Gaussian structure of fractional Brownian motion and its Malliavin derivative process, 
following Section 2 of |l8j. Let S be the set of step-functions on R with values in M*^. 
Consider the Hilbert space TC defined as the closure of £ with respect to the scalar product 
induced by 

d 
i=l 

for any — oo < Sj < < +oo and — oo < ti < f < +oo, and where Rnit, s) is given by 
(!77|) . The mapping 

d 

(l[ti,ti], • • • , l[id,t'']) ^ i^t' ~ ^*«) 

i=l 

can be extended to an isometry between Ti. and the Gaussian space Hi (B) associated with 
B = [B^, . . . , B"^). We denote this isometry hj ip ^ B{ip). Let S be the set of smooth 
cylindrical random variables of the form 

F = f{B{^i), . . . , 5(<^fc)), ^.en, z = 1, . . . , fc, 

where / G C°°(]R'^''^, M) is bounded with bounded derivatives. The derivative operator 
D of a smooth cylindrical random variable of the above form is defined as the H-valued 
random variable 



This operator is closable from L^(f2) into L'p{Q]T-C). As usual, D^'^ denotes the closure of 
the set of smooth random variables with respect to the norm 

WFWl^ = E|F|2 + E||DF||^. 

In particular, if D^F denotes the Malliavin derivative of F G D^'^ with respect to B^, we 
have D^B-l = 5jjl[o,t] for i,j = l,...,d. 

The divergence operator I is the adjoint of the derivative operator. If a random variable 
u G L^(f2;7i) belongs to dom(/), the domain of the divergence operator, then I{u) is 
defined by the duality relationship 

E{FI{u)) = E{DF,u)n, (80) 
for every F G D^'^. Moreover, let us recall two useful properties verified by D and /: 

• If M G dom(J) and F G D^'^ such that Fu G L^(fi; 7i), then we have the following 
integration by parts formula: 

I{Fu) = FI{u) - {DF, u)n. (81) 

• If M verifies EIImH^+EHD-uH^^^ < oo, DrU G dom(/) for all r G M and (/(-Drw))reR 
is an element of L^{Q; Ti), then 

Drl{u) =Ur + I{DrU). (82) 

5.3. Delayed Levy area and fractional Brownian motion. The stochastic integrals 
we shall use in order to define our delayed Levy defined, in a natural way, by 

Russo-Vallois symmetric approximations, that is, for a given process 0: 

1^ d^Bl = - lim {2e)-' 0^ - Bl^^)dw, 

provided the limit exists. This pathwise type notion of integral can then be related to 
some stochastic analysis criterions in the following way (for a proof, see [IJ): 

Theorem 5.1. Fix t > and let G ©^'^(T^) be a process such that 



exists, and such that, setting i{t, u) = u^^^^ + (^ ~ uY^^^ for < u < t, 

i{t, u)dudr < oo. 







/[0,t]2 





Then (l)d°B^ exists, and verifies 

f(f)d°B' = /^'(01[o,t]) + Tr[o,t]D^>. 
Jo 

With these notations in mind, the main result of this section is the following: 

Proposition 5.2. Let B be a d- dimensional fractional Brownian motion and suppose 
H > ^. For V G [— r, 0], let B^{v) be the delayed Levy area given by: 

Bliv) = I dB^® r dBr, I. e. Bl{v){i,3) = [ dB^ [ ^ dB^, t,j G {!,..., 4, 

J s J s+v J s J s+v 



Then almost all sample paths of B satisfy Hypothesis \3.4 



Proof. When H = ^, the desired conclusion is easily obtained, because the Russo-Vallois 
symmetric integral coincides with the Stratonovich integral. Moreover, for H > ^ the 
Russo-Vallois symmetric integral coincides with the Young integral, which is well defined 
in this case, and the assertion still follows easily from the properties of Young integrals. 

Now, fix i < if < i. It is a classical fact that B e C7([0,T]; R'^) for any i < 7 < if. 
Due to the stationarity property (1791) we will work without loss of generality on the interval 
[0, t — s] instead of [s, t] in the sequel. 

1) Case i = j. When t> = 0, it is easily checked that 

E\Blm^,^)\' = \m-B.,\^=^-\t-sr. 

Let us now consider the case where t> < 0. For = {B^^^ — i?*)l[o,t-s](-)) the conditions of 

Theorem 15.11 are easily verified, hence * (j)ud°Bl^ exists. Notice moreover that we have 

Dr'<Pu = '^[v,u+v]{r)l[o^t-s]{u) and, for u E [0,t — s] and e G [0, —v] (which is always the 
case, for a fixed v < and e small enough) it holds 

{ll^^u+v]A[u-e,u+e])H = + £\ " 1^ - ^1 + \v - U - e\ - \v - U + e\ ) 

= ^ {{-V - £)2^ - {-V + ef'' + {-v + u + 5)2^ -{-v + u- ef") . 
Thus, we obtain 

Tr[o,-.]Z^^> = -H{-vr-\t -s) + \{{t-s- " (-^)'^) • 

For X > 0, it is well-known that < ((— f ) + x)^^ — {—vY^ < 2H{—vY^^^x. Applying 
this inequality to the second term of the right hand side of Tr[o,t-s]-D'^V) ^6 get 

Tr[o,_.]D^>| < H{-vf'''\t - s). (83) 
On the other hand, we have by ( !82l) 

Df'J^'(0) = 0, + J^'(Df'0) = (0, + /^'(l[.-.,+oo)n[o,t-.]))l[o,t-.](r). (84) 

When —V >t — s, then [r — v, +oo) fl [0, t — s] = for any r G [0, t — s]. By using 
we deduce 

E\I^ (0)|2 = E||0||^ = E\\B\^ - 5*||^([o_4_^]) = E||fi*||^(p^4_^]) 



where the two last equahties are due to the stationarity fl7^ and scaling flTHl) properties 
of fractional Brownian motion. 
When —V < t — s, then 

I^'il[r-v,+oo)n[o,t-s]) = {Bl_s - 5*_^)l[o,t-s+^](r). (86) 



We deduce 

Ei/^'(0)r 

= E{DI^\<P),^)n by 



— 


E||<^ 


'\\H{[0,t~s]) + (l[r- 


--i;,oo)n[o,t-s]), 0)w([o,i-s]) by 






E\\6 


i||?./r. n + EiiB'' - 


-'-^ — -"-[Oji-s+uJ ; Y^/7i([0,t-sJ) 


bv flHHIl 


< 


E||<^ 


>||^([o,_,]) + E(||(i?,^, 


- -SLi,)l[0,t-s+t,] H([0,t-s]) </ 


^||H([0,t-s])) 


< 


1^1 




-s ~ -^■-D)l[0,t-s+t)] H([0,t-s]) 


because ah < ^(a^ + 6^) 


= 


2^ 




-E||(5j_^+^ - 5*)l[o,t_s+„]|| 


|?,([o,t-.]) by (ITHD and (ITHD 




3. 


- s)'^'^E 5* ^([Q^y) + 




-^(i+s-t))-) H{[0,1]) 




3. 
2** 


- s)'^'^E||5'||^([o^i]) + 


\{t-s + vY"E\\{B\-B^) 


||?,([o,i]) bydZHl) 


< 


-it 

2^ 


- sY" (3E| i?'|^([o^i]) 


+ E||(5i-fiO||^([o,i])). 


(87) 



Finally, we can summarize and fl87p in 

E|/^"(0)P<CH|t-.r, 

with a constant ch > 0, in particular independent of v. Putting together this last estimate 
with inequality (l83l) . we end up with: 

E\BUv)it,t)\'<CHil + \v\'''-')\t-s\''', 

for any v G [— r, 0]. 

Case where i j- By stationarity fl7^ . we have for any v G [— r, 0] that 
f RJ' _ Ri ^ f Ri R« 

Thus, the delayed Levy area Bq j_^(t>)(z, j) = ^ {bI^^ — Bi)d° Bl^ for t> < behaves as in 
the case where v = 0. But it is a classical result that Bq j_^(0) is well-defined for H > 1/3 
(see, e.g., [20]). Moreover, it follows again by the stationarity (I79|) and the scahng ( ffSl) 
properties that 

EiBoV.(^)(^, j)r = E|BoV.(o)(^, j)r = \t - sr^E|Bo^(o)(., j)r < c^.it - sr. 

Immediately, we deduce that 

E\BUv){^,J)\'<CH\t-sr 

for any v G [— r, 0]. 

Both in the cases i = j and i ^ j, the substitution formula for Russo-Vallois integrals 
easily yields that 5B^(f) = SB"" ® 55. Furthermore, since B^(w) is a process belonging 
to the second chaos of the fractional Brownian motion B, on which all norms are 
equivalent for p > 1, we get that 

E\Bliv)i^, j)f<c,\t-sf^'' 



for i ^ 2 ci-iid 

E|/^'(0)|P < Cp|t-spP^ (89) 

when i = j. In order to conclude that B^(f) G €2'^ {M.'^^'^) for any | < 7 < and 
V e [— r, 0), let us recall the following inequality from pj: let g E C2{V) for a given 
Banach space V; then, for any k > and p > 1 we have 

\\gh<c{U,+2/pA9) + \\Sgh) with UrA9)=(l I Tf^dsdt] . (90) 



By plugging inequality dHED-dHD into by recalling that 5B'^{y) = ®5B and 
hold, we obtain that B'^{v){i,i) e Cl^ iW^""^) for any i < 7 < and j = 1, . . . , c?. □ 
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